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^ ■ Abstract 

0^ ; 

^ ■ Quasielastic antineutrino/proton and neutrino/neutron scatterings can be well ap- 

^ . proximated by simple formula3, valid around MeV or GeV energies. We obtain a 

I single expression valid in the whole range, and discuss its relevance for studies of 

■ supernova neutrinos, which reach intermediate energies. 

O 

(N ; 

cn ! 1 Introduction 

O ; 

r~| ! The quasielastic reaction v'ep — > e~^n (inverse beta decay) is of special relevance for studies of Ue at 
Qh' sub-GeV energies P]. Indeed, (1) its cross section is relatively large^, ~ GpE^; (2) it can be accurately 
Q . computed; (3) it has a low threshold > 1.806 MeV rUe + mn — nip; (4) the measurable e"*" energy 
|H \ is strongly correlated with the i>e energy; (5) materials rich in free protons are cheap (e.g. water, 

■ hydrocarbon) and this permits to build large detectors; (6) in scintillators, it is possible to tag both 
. . \ the e"*" and the neutron, reducing the backgrounds. These features are not shared by other reactions, 

' 6.g. (1), (4) and (6) do not hold for ve ve (since mg <C m^); Ve^^O — > e''~^^N has high threshold 
^ . Ey > 11.4 MeV and it is theoretically much less clean. 

^ \ The cross sections for quasielastic neutrino and antineutrino scattering on neutron or proton have 

■ - - ' been calculated long ago for applications at high energy Ey ~ GeV (e.g. accelerators). The reference 
formulae for these applications are eqs. (3.18) and (3.22) in the review of Llewellyn-Smith routinely 
used, for instance, in analyses of atmospheric neutrinos. Vogel and Beacom ^ remarked that effects 
suppressed by (m„ — mp)/Ey have been neglected in eq. (3.22), which is a bad approximation at low 
energy Ey ~ MeV (e.g. z/g generated by nuclear reactors or supernovae). These authors proposed 
theirs eqs. (12) and (13) which, neglecting effects suppressed by powers of Ey/rup, provide a simple 
and precise description of neutrino/nucleon scattering at low energy. In this paper, we derive a unique 
formula valid at low, intermediate and high enGrgics, without assuming ttlyi — *^ Ey nor Ey <^ nip. 
We will compare our cross section with other approximations of common use, and discuss the range of 
energies where they are reliable. Our findings are relevant in particular for (1) detection of supernova 
i>e, that can have an energy between few MeV up to about 60 MeV; (2) transport of e^, i^e and i>e 
inside the collapsing star, where the effective temperatures are higher and neutrino energies can reach 
200 MeV. 

The cross-section is computed in section |^ and discussed in section |31 

^It is interesting to recall that the first authors who computed it concluded; ' This meant that one obviously would 
never be able to see a neutrino. ' 
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2 Cross section 



We write explicit expressions for inverse beta decay (IBD), namely quasielastic scattering of electron 
antineutrino on proton: 

I'eiPu) +p{Pp) e+(pe) + n{pn) IBD reaction. (1) 

We indicate how to obtain the corresponding expressions for other antineutrinos and i^r, and for 
quasielastic neutrino/neutron scattering. Let us define 

A = m„ - mp « 1.293 MeV, M = (m^ + m„)/2 938.9 MeV. (2) 

The differential cross section at tree level in the weak interactions, averaged (summed) over initial 
(final) polarizations is: 

da _ G^cos% ,^^2 1 .o^ 

M ~ 27r(s-m2)2 I ^"^^ 

where Gp = 1.16637 10~^/GeV^ is the Fermi coupling extracted from /i-decay, cos^c = 0.9746 it 
0.0008 U] is the cosine of the Cabibbo angle and A4 has the well-known current-current structure: 

M = Uve7"(l - l5)Ve • Un (^flja + 5l7a75 + «/20-„6^ + ^2-^75^ Up (4) 

where q = Pu — Pe = Pn — Pp- A straightforward calculation yields^: 

\M^\ = A{t) - (s - u)B{t) + {s- ufC{t) (5) 
where s = {py + Pp)'^, t = [p^, — Pe)'^, u = {py — pn)^ are the usual Mandelstam variables and 



16 ^ = (t - ml) 



-A2 



4|/i2|(4M2 + t + ml) + A\gf\{-AM^ + t + ml) + \f^\{t^/M^ + 4t + Aml)+ 

+Amlt\gl\/M^ + 8Re[/i*/2](2t + ml) + lQmlRe[glg2] 

(4|/2| + t|/||/M2)(4M2 +t-ml) + 4|<7f |(4M2 - t + ml) + Aml\gl\{t - mD/M^ (q) 

+8Re[/i*/2](2t - ml) + lQmlM9l92] - 32m2MARe[gt (A + /s)] 
IQB = 16tRe[<7i*(/i + /a)] + 4m2A (|/|| + Re[/f/2 + 2glg2]) /M 
lQC = ml\ + \9l\)-t\f^\/M\ 

In the limit A = eq. (jUJ reduces to eq. (3.22) of Llewellyn-Smith. 

Since strong interactions are invariant under T and C, the adimensional form factors fi,gi are 
real functions of the transferred 4-momentum t = < 0. As discussed in section 13.21 they can be 
approximated with |5] 

_ {l-(l + g)t/4M^e} gi(0) _ 2M^gi 
^^^'^'^ (l-t/4M2)(l-t/M^)2' ^1 (l-t/M2)2' ^2 ^2_i 



^We give some details of the computation. The polynomial structure in s — it can be easily understood: if we express 
the leptonic tensor L""^ in terms of Pi — Pv + Pe and g, and the hadronic one in terms of Ph = Pp + Pn and q, it is 
sufficient to note that pi ■ ph = s — u (while = t, q ■ pe = —ml and q ■ Ph = rnj^i — rn^ = 2MA). The dependence on A 
of \M'^\ can be understood by deepening the same argument, namely: (1) one decomposes the leptonic tensor L"'' in 4 
parts: p^p^, q'^q^ — q^g"^^ , rn^g"''', e'^^'^'^pia<id\ (2) then, one considers the general Lorentz-invariant decomposition of the 
hadronic tensor H{ph, q) (6 parts); (3) finally, one checks explicitly the 24 contractions, and verifies that A appears only 
in the combinations A'^ and m?A. 
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where 51 (0) = -1.270 ± 0.003 0, = 0.71 GeV^, « 1 GeV^. Finally ^ = Hp - Kn = 3.706 is 
the difference between the proton and neutron anomalous magnetic moments in units of the nuclear 
magneton (k^ = 1.792, k„ = —1.913). 

The analogous expressions for ly^n scattering are obtained by exchanging pi, <^ —Pe, and nip — > m„ 
in the flux factor. Therefore, \A4'^\ keeps the same expression with s — u ^ u — s; in view of its explicit 
form, this amounts to just replace B — > —B in eq. ©. The expressions for or Vr scatterings are 
obtained by replacing me with or m,-. 

Our final expression is somewhat involved but allows to numerically compute any desired integrated 
cross-section. It is also useful to present some simple analytic approximation that is valid at low 
energies (say, at supernova and reactor energies or below). In fact, the fg, z^e cross sections simplify, 
when expanding in powers of the parameter: 

e = Ejmp (8) 

When we prescribe that E'g, A, rrif, ~ e, we get 

s — rrip'-^e, s — u^e, t . 

The flux factor in the denominator of da/dt gives a 1/e^, and the amplitude |A^^| starts at order . 
The leading order (LO) approximation in e is obtained by retaining only 

A ^ M^ifl - gl){t - ml) - M^/\\fl + gl) 

B ~ (9) 

Since i? ~ at LO t'gn and v^p scatterings are described by the same Ai^. 

A simple approximation which is accurate enough to describe the detection of supernova neutrinos 
events via IBD is obtained using the expression of M."^ accurate up to NLO in e, given by 

A M\fl - gl){t - ml) - M^A^fl + gf) - 2mlMAgi{h + /a) 

B^tgi{fi + f2) (10) 

C^{fl + 9l)/4. 

The t dependence of the form factors contributes at next order (NNLO): at NLO the form factors are 
approximated with constants. 

2.1 Total cross section, da/dE^ and event numbers 

The cross section expressed in terms of the neutrino and electron energy in the rest frame of the 
proton. El, and E^, is particularly useful. Inserting 

s — mp = 2mpEy, s — u = 2mp{Ey + E^) — ml, t = — — 2mp{Ey — E^) 

it is given by 



dhe dt Im. 



-{E,, Ee) = 2mp— if E, > ^thr = ^ ^ (11) 

The allowed values of E^,, Ei < E^, < E2, correspond to the possible scattering angles 6'^^ in the 
center of mass (CM) frame: 

2 2 2 

Ei,2 = E,-5- —E^'^iEr ± pr), with 6 = UhLZ^hlZ^ (12) 
mp 2mp 
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Figure 1: Fig. QJa.- total cross sections of quasielastic scatterings at supernova energies. Fig. Q).- 
average cosine of the charged lepton scattering angle (cos 6) . 



where the energy and momenta in the CM have the standard expressions: 

,2 



mt 



2^s ' 



rni + -mi 



2^ 



(mn+me 



2^ 



(13) 



For the neutrino reaction, one has just to replace run rnp in all formulae above except in |A^^|, 
where A remains positive, with the proviso to set E'thr = for (it would come negative). The total 
cross section is plotted in figure^ and tabulated in tabled We have included also two corrections: 



1. the radiative corrections to the cross section [HIIZI, well approximated as [Zj 

3 



(14) 



where a is the fine-structure constant. 

This one-loop correction is valid at E^ <C mp, when the final state positron (electron) is the 
only radiator, and in the hypothesis that all the energy Ef. in electrons and bremsstrahlung 
photons is detected. In general the differential cross section is given by a more complicated 
expression which depends on detection threshold for photons. The total cross section is still 
given by eq. ()14() . if all final state electrons are detected j7j. 

High-energy electroweak corrections are automatically included in the values of 51(0) and Gp 
extracted from low-energy experiments. 



2. We multiply a{y^n pe) by the Sommerfeld factor that accounts for final-state interactions: 



F{Ee) 



1 — exp(— 7?) 



with 7] 



(15) 
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Table 1: Values of a{V(,p — > ne"^) in 10 cm^ and of the average positron energy and cosine of the 
scattering angle. All energies are in MeV. 

Both the radiative corrections and final state interactions amount to a ~ 2% correction, which can be 
important for high statistic data samples. 

The formulae given above can be adapted to the case when the nucleon is not at rest (e.g. inside a 
supernova). Similarly, one can include Fermi blocking- factors in final state to account for degenerate 
electrons sea in the supernova. 

It is straightforward to compute the average values (x) = J x da/ J da of interesting kinematical 
quantities x, like x = E^, El — {E^)'^, cosO, etc. We plot in figure 2 and tabulate in tabled the 
average cosine of the scattering angle of the charged lepton. Detecting supernova neutrinos via the 
IBD reaction, one could try to track their source by measuring the direction of positrons, that are 
copiously produced |H1|3] (more discussion below). Another important quantity is the average lepton 
energy, which can be approximated (at better than 1% below 100 MeV) by 

(Ee) « (El + E2)/2 = E,-5- Ef'Ef'lm^ (16) 

This is a much better approximation than the usual formula Ef, = Ey — A; it permits to relate E^. 
with Ey incorporating a large part of the effect due to the recoil of the nucleon. 

It is straightforward to include cuts or efficiencies: e.g., if the positrons are detected only in the 
window El < Ef. < Eh we define 

a{E,) = / dEe -pEriE,, E^) ^{E, - Ei)^{Eh - E^) (17) 

J El drje 

where 'd{x) is the step function, equal to (1) for x < {x > 0). The total number of events in the 
detector and the distribution in positron energy are given by 



f ,^ d^ , dN„+ f d^ da 



5 



where Np is the number of target protons and $ is the z/e flux. We do not have a simple formula for 
the exact limits of integration in dEy. A simple recipe is to integrate numerically in the larger range 

Ee + 5<E,< ^ft^J,,/ (19) 
1 — 2{Ee + d)/mp 

obtained replacing and E^^ with E^^ in eq. ()12|1. enforcing the correct kinematical range by 
means of ??-functions; the upper limit has to be replaced with infinity when the denominator becomes 
negative. 



2.2 Angular distribution 

The cross section differential in the angle 9 between the neutrino and the positron is obtained from 
da/dt after including the Jacobian coming from t = — 2Ey{Ef. — cos 9): 

d(J , ^ Ve£ da , , 

iE,.cos9) = , , -=r, (20) 



dcos9' ' i + e(i_|^cose) dE^' 

where e = Eyjmy. In eq. (|2U() . E^, and Pe are functions of and Q'. 

j,^=iE^-m^e)^ecos9^[K-^Y-mi.^ p. = T^T^ (21) 



with K = (1 + e)^ — (ecos^)^. A lepton of low energy cannot be emitted in all directions 9 = O-^vr, since 
it is "dragged" by the motion of the CM.^ The r, being heavier than the proton, is always emitted in 
the forward cone. For z/g or a non trivial allowed range of 9 is obtained only in a very small region 
energy range just above the threshold -Ethr (given in eq. (jllf) ): 

, [nin — [nip — m-e)]^ 

Ethr < Ey < E^^ = £/thr + 'me — 7 ^ — (22) 

2mp[mp — me) 

For the muon (electron), E[^^^ — i?thr = 0.8 MeV (1 eV). Above E^^^. all scattering angles are allowed. 
Therefore, we can conveniently simplify I'f., i^^, scatterings by replacing E'thr with E'^^^j. and considering 
the full range 9 = ^ tt. This is a very good approximation. 



3 Discussion and applications 

How important is to use an accurate cross section, or in other words, what is the error that we 
introduce when we use an approximate cross section, instead of the correct one? The answer, of 
course, depends on the process we are considering. Here we compare our cross section with other ones 
that can be found in the literature, evaluate its uncertainty, and finally discuss the specific application 
for supernova signal. 

Here we mention some other possible applications of these results. A precise cross section can be 
used in an accurate description of neutrino transport in a collapsing supernova (see e.g. [HI for some 

^The region of "CM-dragging" is defined by the condition that the CM velocity /3 is larger than /?^". In this region, 
the function 9'^" (9) — or 6l°"(f) or d^'^{Ee) — is not single- valued: there are two values of 9^'^ that lead to the same 9. 
This is evident from the relation tan 61 = sin 6l*="/[7(cos ei^" + /^//J^")] where 7 = l/y'l - (3'^. The region of 'dragging' 
is visible in the leftmost part of figure 2: it is where the curve "i/e on proton" rises suddenly. Exactly at the threshold, 
the positron is emitted only at 9 = 0, so that {cos 9) —^ 1. 
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Table 2: Percentage difference between our full result and various approximations for v^. (above) and 
i^e (below) total cross sections. A negative (positive) sign means that a certain cross section is an 
over(under)-estimate. It is easy to implement approximations maked with -kirir, while implementing 
those marked with a -k is not much simpler than performing a full computation. 

recent calculations). We recall that neutrino heating (or more generally, the coupling of neutrinos 
with the deepest layers of the mantle) takes place at energies E^ of several tenths of MeV, and it is 
regulated by the cross sections that we are discussing. Furthermore, electron neutrinos with energies 
in the range 50 — 150 MeV are particularly important for the dynamics of neutrino emission: they 
are around the Fermi surface of the degenerate neutrino distribution, and determine the neutrino flux 
from the inner core. 

A precise cross section should also be used to interpret searches for U(. or i/^ low energy neutrinos 
from the sun (see e.g. ^Oj), from gamma-ray-burst or other presently unknown cosmic sources or 
also, in studies of z/^ or interactions just above the threshold of /x production (e.g., for large sample 
of data in terrestrial experiments). 

The LSND and Karmen experiments |12j used a flux with a continous spectrum up to 52.8 MeV, 
mainly produced by decay at rest, looking for appearance using the i>eP e^n reaction. The 
experimental result is controversial. If such a signal exists, a precise computation of the VeP cross- 
section would be necessary to interpret a more precise future experiment. 

3.1 Comparison with various simple approximations 

Table m shows the percentage difference 

5a = 100 X (1 - a' /a) (23) 

between our cross section a and various approximations a' , evaluated at the same input parameters 
(percentages below 0.1% have been set to zero). The first row shows the neutrino energies at which 
the comparison is made. In the upper panel we consider the reaction; in the lower panel the i^e 
reaction. 

1. The naive low-energy approximation (see e.g. jl3j ) 

a w 9.52 X 10"^^^^^ cm^ Ee = E^ ± for Ve and v^, (24) 



7 



obtained by normalizing the leading-order result to the neutron lifetime, overestimates a{i>eP) 
especially at high energy. It is not recommended for analyses of supernova neutrinos, nor for 
precise studies of reactor Uf,. The naive a{i'^n), instead, agrees well with the exact cross section. 

2. A simple approximation which agrees with our full result within few per-mille for Ei, < 300 MeV 
is 

a(PeP) ~ 10-^3 ^^2 p^^^ ^-0.07056+0.02018 In £.-0.001953 ln3 ^ = - A (25) 

where all energies are expressed in MeV. 

3. The low-energy approximation of Vogel and Beacom jlj (which include first order corrections 
in e = Ey/nip, given only for antineutrinos) is very accurate at low energies (£^,^^60 MeV), 
however underestimates the number of supernova IBD neutrino events at highest energies by 
10%. Higher order terms in e happen to be dominant already at £'jy>135MeV, where the 
expansion breaks down giving a negative cross-section |^ . 

4. The NLO low-energy approximation, defined by eq. (jlUj) . can be used from low energies up 
to the energies relevant for supernova i>e detection. (We expand the squared amplitude in e 
but, unlike Vogel and Beacom we treat kinematics exactly, so that some higher order terms are 
included in our NLO cross section). 

5. The high-energy approximation of Horowitz (l^j . obtained from the Llewellyn-Smith formute 
setting vrie = 0, was not tailored to be used below 10 MeV, and it is not precise in the 
region relevant for supernova neutrino detection; however, it is presumably adequate to describe 
supernova neutrino transport. 

6. The Llewellyn-Smith high-energy approximation, improved adding m„ 7^ mp in s, t, u, but not 
in is very accurate at all energies relevant for supernova neutrinos, failing only at the lowest 
energies. As proved in section |21 this is a consequence of the absence in of corrections of 
order A/mp. 

7. Finally, approximation 6. can be improved by including also the dominant low-energy effects in 
the amplitude Al, as discussed in section IIB of With respect to our full result, this amounts 
to neglect A and mg in A, B, C, eq. ©, and reinsert part of them modifying A — > A — AM^A'^C. 
This approximation agrees with our full result within few per mille. 

3.2 Overall uncertainty 

We now discuss how accurate our full expressions for the cross sections are. All input quantities have 
been precisely measured, except the t-dependence of the form factors of eq. ((TJ. As shown in section 2, 
it affects the cross sections only at NNLO in Ey/nip. In practice, taking the form factors as constants 
overestimates the total 9f.p and cross sections by about 1% • (£'i,/40MeV)^. At low energy this 
is an excellent approximation, and the dominant uncertainty is due to Oq and 51(0). The values we 
adopted are extracted from different sources. The Cabibbo angle Oq is extracted independently from 
Vus (measured from i^£3 decays) and from Vud (measured from nuclear decays dominated by vectorial 
matrix elements, and from neutron decay). Since the two determinations agree only within 2.2a and 
theoretical errors play an important role, the total error has been conservatively increased [5]. 

The axial coupling (71 (0) is measured from neutron decay^. Different experimental determinations 
do not fully agree, therefore we conservatively increased the error. Newer measurements, performed 

*Its uncertainty is fully correlated with the uncertainty on the 9c determination from neutron decay, since neutron 
decay measures the same combination of (?i(0) and 6c that enters in IBD. 
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with a higher neutron polarization than older ones, are consistent and agree on (7i(0)//i(0) = —1.272 it 
0.002 when older determinations are discarded — a value slightly different from the one quoted in 
section 2. Isospin-breaking corrections to /i(0) = 1 are negligible |15j . 

In conclusion, at low energy a{i'ep) has an overall 0.4% uncertainty, which is adequate for present 
experiments. The ratio between the measured and the no-oscillation reactor v^, flux is 1.01ib2.8%(stat)ib 
2.7%(syst) at CHOOZ |16j . The KamLAND collaboration presently has a 6% systematic error on their 
reactor data ^7]i ^"^^ aims at reaching the 4% level in future. Concerning supernova z/g, experi- 
ments should collect few hundreds to several thousand of events, depending on when and where the 
next galactic supernova explosion will occur pi 1181 IT^ . 

We finally discuss the uncertainty introduced by the t dependence of the form factors, which for 
the largest supernova i>e or i^e energies is not fully negligible. The SU(2) isospin symmetry (broken 
only by small quark masses and by electromagnetism) allows to relate the vector form factors /i and 
/2 to measured eA^ cross sections, and conservation of the vector current predicts /i(0) = 1. The first 
terms in the expansion in t 

flit) ~ 1 + 2M/ GeV2, /2(t) ~ e[l + 3.4t/ GeV^] 

are known with about 1% accuracy [20j . The phenomenological parametrization in eq. Q, partially 
suggested by theoretical considerations, reproduces reasonably well the experimental data. Employing 
slightly more accurate parametrizations for the vector form factors /i and /2 (e.g. j20j l would not 
improve the accuracy of our results due to the larger uncertainty on the axial form factors gi and §2- 

The axial form factor g2 (written in eq. ((TJ in terms of gi , as suggested partial conservation of the 
axial current) has a negligible impact on i^g and i/g scattering at intermediate energies. Thus, the main 
uncertainty comes from gi. Vector meson dominance would suggest gi = gi{0)/{l — t/M\^) where 
M^j = (1.23 ±0.04) GeV is the mass of the pseudo- vector that, mixing with the W boson, should give 
the dominant effect. The parametrization of gi adopted in eq. ((T)) is obtained by squaring (1 — t/M^) 
and leaving Ma as a free parameter because this provides a better fit of data. Lattice simulations do 
not contradict this hypothesis within the present theoretical errors, that are however still large |21j . 
Assuming that the parametrization is correct, one obtains 

1. Ma = (1.069 ±0.016) GeV |12 from vr electroproduction on nucleons (7*A^ ttN'). The 
presence of tt makes this probe experimentally precise but theoretically uncertain, as we are 
interested in {n\ - ■ ■ \p), not in (7rn| ■ ■ ■ \p). 

2. Ma = (1.026 ± 0.021) GeV from global fits of z^A^ scattering data. Most of the data are 
not taken at |t| <C GeV^, where the parametrization of gi is certainly justified. 

3. /i capture on p could provide a clean and direct probe of g2 at t = — O.SSm^. The present 
experimental data |24j hint to a value about 50% larger than the one we assume^, but their 
interpretation is under debate • New data should clarify the issue j25j . 

The above discussion shows why it is difficult to assess the uncertainty on gi and g2- Optimistically 
assuming that 1. or 2. is right, it is negligible. On the other side, a pessimistic estimate can be 
obtained by using Ma-^ in place of Ma- the total I'eP cross section increases by 0.4%(£'j,/50 MeV)^ for 
El, < 200MeV. The shift remains relatively small because, as shown in section 2, the t-dependence of 
the form factors affects UeP only at NNLO in Ey/rup. 

^Such an increase of g2 has some impact on r/^ and scattering, and would decrease a{u^p) by 3.8% x (150 MeY / E,y)^ . 
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3.3 Search for relic supernova neutrinos 



Here we discuss the apphcation of our cross section to the search for neutrino radiation from past 
core-cohapse supernovae. We focus on the energy range Eg > 20 MeV which is particularly suited 
for this type of search, since it contains a large part of the expected signal and it is largely free from 
background (for an experimental proof, see EHl)- We assume that the energy distribution of supernova 
in the detectors is described by the following differential flux (see e.g. |Ul 

d$ Etot 



dE„ 47rZ)2 



cos 6112 ni^— j + ^ sm Ou n-l^— 



(26) 



where 1) 612 ~ 34° describes 'solar' oscillations in supernova mantle; 2) n{x) = 120x2/77r^[l + e^] is 
the Fermi-Dirac occupation factor; 3) fe ~ 1/5 and /^^r ~ 1/7 are the energy shares in i/g and in 
i^fj,,T] 4) Te ^ few MeV is the temperature of Uf, (before oscillations) and 5) T^j^^j- is the temperature 
of I'fi^T, tliat we assume to be T^^,- ^ 1.3 x Tg; 6) E'tot is the total energy released in neutrinos, D 
the distance of the supernova. We neglect various effects that would not change our conclusions, 
like 'pinching' of the spectrum (i.e. non-thermal effects), eventual oscillations in the Earth, existence 
of other detection reactions (we focus on detectors with a large fraction of free protons like water 
Cerenkov or scintillators), precise definition of the search window (which depends on the background 
level), detector efficiencies, other oscillations eventually induced by a non- vanishing ^13, etc. 

We compare the expected number of events A'^g+ in the energy range E^ = [20, 40] MeV using (a) 
the cross section of eq. H17|) and (b) the naive (and much used) cross section in eq. (|24() . The latter 
leads to an overestimate: 

/ N,+ (h)\ ^ ^ , Te- 4.5 MeV , , 

100 X 1 - J ) ' w -25% + 2.4% X X, where x = " . 27 

V N^+{a)J MeV ^ ' 

The Super-Kamiokande collaboration presently gives the dominant bound on the relic Ve flux, con- 
strained to be less than 1.2i>e/cm2s at 90% CL for E^ > 19.3 MeV [201 — slightly ab ove the flux 
expected by cosmological models |27|. Since this bound has been obtained assuming the naive u^P 
cross section, eq. (P7|) suggests that it could be weakened by ~ 20% (a reanalysis of the data with a 
precise cross section would allow to be more precise). Similar considerations apply to certain calcu- 
lations of the expected event number from next galactic supernova (e.g. ref. ^Hl) and to estimates of 
the total energy Etot from SN1987A neutrino data (e.g. ref. [lHl)i although the low energy cuts are 
typically smaller than 20 MeV and therefore the inaccuracy is less severe. 



3.4 Low energy atmospheric neutrinos as a background 

Low energy atmospheric neutrinos constitute a background for searches of relic supernova neutrinos. 
At Super-Kamiokande this background is mostly due to muons produced below the Cerenkov 
threshold in water (S^ < 160 MeV), whose decay produces a positron. Assuming that at low energy 
the v'^ flux scales as ^ l/E^, , the cross section with nip = rrin = M overestimates the background 
only by 1.5%. This is smaller than what suggested by a simple minded estimate, performed looking 
at the total cross section without taking into account that the condition that the muon is produced 
below the Cerenkov threshold does depend on A = m„ — rrip. Similar considerations apply to the De 
atmospheric neutrino background: e.g., the use of the Llewellyn-Smith cross section for Ef, = [50,80] 
MeV introduces a 2.4% overestimate. 



3.5 Distributions of positrons from supernova neutrinos 

Now we compute the energy and scattering-angle distribution of positrons generated by supernova i>e 
and detected via the IBD cross section, that we compute according to the formulae of sections 12. H and 
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0.7 




Figure 2: Fig.^jp, shows the energy spectrum of positrons. Fig.^j) shows the positron scattering angle 
distribution, for different lower cuts on the positron energy. We assume that supernova have the 
flux in eq. \2f^) and are detected via the UeP — > ne'^ reaction. 

12.21 Again, we assume the Ue flux of eq. (|26)) . where we set = 4.5 MeV. The result is shown in 
fig. [5J The average cosine of the scattering angle asymmetry is approximatively given by 

(cos 6*) = 0.024 + 0.013 x (28) 

(x as in eq. (|27|) ) which means that the scattering is slightly forward. The forward/backward asym- 
metry is not particularly large, but could be observed in a large sample of neutrinos from a galactic 
core-collapse supernova using wisely designed selection criteria (such as an energy dependent angular 
analysis). For example, it is sufficient to increase the lower cut on the positron energy to increase the 
asymmetry considerably 

r 0.028 for Ee > 10 MeV 
(cos 6') = < 0.055 for > 20 MeV (Te = 4.5 MeV). (29) 
[ 0.102 for Ee > 30 MeV 

as more precisely illustrated in fig. [^d. A lower cut on Eg increases (cos 6) for two different reasons: 
because the IBD reaction is forward peaked at higher energy, and because forward positrons are 
more energetic. (In the Super-Kamiokande detector a lower cut E^'^S MeV must anyhow be imposed 
in order to cut low energy backgrounds). Two remarks are in order: 1) An accurate description 
of the direction of IBD events is also important to separate other classes of events, as those due 
to reactions with Oxygen nuclei in water Cerenkov detectors. 2) An accurate description of the 
spectrum is particularly important when trying to extract information from the high energy tail of the 
distribution, about e.g. oscillations in the Earth, or non-thermal corrections to the energy spectrum. 

Of course, not only (cos 0) but also the number of events A''g+ depends on the temperature. As- 
suming unit efficiency above Eg = 8 MeV and using eq. (|18|) . we find that varying Tg around the 
reference temperature, the number of events varies approximately as l + 0.24x — 0.02a;^: Nf,+ increases 
(decreases) by 23% (26%) for Tg = (4.5 ± 1) MeV. 
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